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ABSTRACT 


The operation of replenishment at sea is investigated for three 
supply ships and L combatants using queueing theory concepts and a 
random walk model in three dimensions. The distribution for total 
replenishment time given the initial number of cambatants to be 
replenished by each supply ship and a specified cyclic order of 
replenishment is expressed in terms of its Laplace transform under 
the assumption of independent, exponential service times for each 
supply vessel. A method for counting all possible sequences of 
replenishment is not found, but same preliminary counting techniques 


are developed which may be useful in its eventual determination. 
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I. INTRODUCTION 


A major aspect of modern naval policy is the stationing of 
task forces throughout the world, and the supply and replenishment 
of these task forces is a critical problem with which naval plan- 
ners are faced. It is vital for ships at sea to be provided with 
a reliable flow of supplies for them to maintain a continuous 
readiness posture. 

Forces afloat are presently supplied with fuel, food, ammunition, 
and other bulk goods primarily by large supply ships which accompany 
task forces during their operations and may periodically proceed to 
bases ashore to replenish themselves. Transfer of goods at sea is 
largely accomplished by cambatants coming alongside these support 
ships and receiving supplies by use of various rigs. The nature of 
the primary alongside method of transfer and its consequent grouping 
of numbers of ships over a period of time increases the vulnerability 
of all participating vessels. Therefore, time of replenishment of 
the whole combatant task force is of primary concern. The problem 
of scheduling a replenishment operation of a combatant task force 
by a group of supply ships such as oilers, ammunition ships, and 
other support vessels has been studied by McCullough [3], Gordon 
and Copes [2], Patterson [5], and Waggoner [6]. Waggoner built a 
mathematical model to study the scheduling of a fixed number of 
combatants by two supply ships and computed the distribution of total 
replenishment time of all combatants. It is the purpose of this paper 


to consider the analogous problem with three supply ships. 


II. FORMULATION OF THE PROBLEM 


In general the problem may be considered fram a queueing theory 
approach with each supply ship representing a service facility and 
each combatant a customer. It may or may not be necessary for each 
custamer to visit each facility, depending upon the needs of the 
custamer and the number of duplicate facilities. Waggoner [6] and 
Milch and Waggoner [4] employ a random walk model for study of a 
two facility queue which operates in parallel and series simultaneously 
with a finite number of customers. Such a model is applicable to the 
replenishment situation. It is the purpose of this paper to extend 
the random walk model to the case of three supply ships under the 
following assumptions: 
(i) A finite number of combatants are initially present, 
and each must be served once by all three supply ships. 

(ii) The service times at each supply ship are independent, 
exponential random variables with parameters \, u, and 
v for ships I, II, and III, respectively, and each one 
serves only one combatant at a time. 

(111) Service is simultaneous and consecutive. When a combatant 
has been replenished in one queue he moves to the next 
queue in cyclic order. Each facility continues in operation 
until it has served all cambatants. 

(iv) No additional combatants join the queves from outside 

the system after the canmencement of operations. 
(v) When a cambatant has been served by all three supply 


ships he departs the system. 


(vi) Transit times between queues are considered negligible. 
Figure 1 illustrates the initial positions of all vessels and 


the order of flow from one queue to another. 
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Figure 1] 


Tit. THE THREE DIMENSIONAL MODEL 


Milch and Waggoner [4] have represented an analagous two facility 
queueing system by a randam walk in a two dimensional lattice. The 
randam walk begins at the origin at time zero and proceeds by horizontal 
and vertical steps to the point (L,L) where L=M+N, and M ard N are 
the initial number of custamers in queues I and II, respectively. A 
horizontal step corresponds to completion of servicing in one queue, 
and a vertical step is a campleted service in the other. Howevera 
all paths to (L,L) are not possible since if X(t) and Y(t) represent 
the amount of custamers served in each respective queue at time t, 
then X(t) < M+ Y(t) and Y(t) < N+ X(t). Figure 2 shows a path in 


such a representation. 


X(t) = L 





M 


Figure 2 


In the two dimensional model the distribution of total operation 


time, T, 1s found in terms of its Laplace transform. 
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If three supply ships are present as in Figure l a three 
dimensional random walk may be formulated which begins at the 
origin at time zero and proceeds by unit steps at a time in one 
of the three axial directions to the point (L,L,L) where L=M+N+P, 
and M, N, P are the initial number of ships in the first, second, and 
third queue, respectively. All paths to (L,L,L) do not represent a 
replenishment operation, however, since if X(t), Y(t), and Z(t) are 
the number of ships served in queues I, II, and III respectively, 
then X(t) < M+ Z(t), Y(t) < N+ X(t), and Z(t) < P + Y(t) are 
constraints that must hold during the whole operation. These con- 
straints define boundary planes for the three dimensional lattice. 

The random walk paths are further bounded by X(t) =L, Y(t) = L, 
Z(t) = L. For the purpose of clarity in further development these 


planes and their intersections will be catagorized as follows: 


Boundary Constituents 
A X=2Zt+M, z=L, and their intersection 
B y=xt+N, y=L, and their intersection 
G Z=y+P, z=UL, and their intersection 
D Intersection of A and B 
E Intersection of A and C 
F Intersection of B and C 
Note: D is not regarded as part of Aor B. Similar statements may 


be made about E and F. 


ime 


IV. DISTRIBUTION OF TOTAL REPLENISHMENT TIME 


A. PATH PROBABILITIES 

The first move of the random walk is made one unit in an axial 
direction with probabilities AJA tuty, ufrituty, and JA +tuty 
for axes X, Y, and Z, respectively, since these are the probabilities 
of stations I, II, or III, respectively, completing replenishment of 
the first combatant under the assumption of independent, exponential 
service times. Due to the memoriless property of the exponential 
distribution these probabilities are true at all steps along the randam 


walk, except at points on the boundaries whose probabilities are as 


follows: 
Boundary Axial Direction Probability 

A ng u/u tv 
Z v/u + v 

B xX AJA *V 
Z V/A + v 

e X vat i 
AG u/rA + yu 

D Z one 

E Y one 

re xX one 


The total time for the replenishment is the time needed for the 
random walk to proceed from the origin to the point (L,L,L) which is 
the sum of the times needed to make each of the 3 L steps required for 


the walk. The distribution of time for any of the steps not originating 


nz 


on a boundary is the minimum of three independent, exponential 
distributions and therefore is itself exponential with parameter 
h+wp+v. Similarly the duration of steps originating on bound- 
aries A, B, and C is exponential with parameters un + v , X} + v, and 
Atu, respectively. The distribution of steps originating on bound- 
aries D, E, and F is exponential with parameters v, u, and i 
respectively. 

Let Q(i,j,k,%,m,n) be a route that has exactly 1, j, k, 2, m,n 
positions on boundaries A, B, C, D, E, F, respectively. Then 
Q(i,j,k,2,m,n) has 3 L - (i+ j +k +2+m+n) positions along the 
path that are not on a boundary. The total replenishment time, T, 
is composed of the following random variables, given that Q(1,j,k,2,m,n) 


@6Curs ; 
ge Be. ae Gy, ee (4-1) 


where A. is the sum of i steps each of which originates on boundary A, 
~stae G Ce Do. EY and FL are Similarly defined. oy is the sum of 

g steps none of which originates on any boundary where g=3 L-i- 

7 —-& - 2 -m -wM. 

Path Q(i,j,k,2,m,n) contains L steps in each of the three axial 
directions. If the L steps in the x direction are examined, exactly 
n of them will begin on boundary F each with probability one. Some 
number, Wye will begin on boundary B each with probability \/d + v 
and some number, k,, begin on boundary C each with probability A/d + u. 
The remaining L-n- j ikem ky steps will begin from a non-boundary 
position each with probability A/A} +u+ v. Similar statements may be 
made concerning the steps in the other two axial directions, and the 


results are tabulated below. 
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Axis Number Step 
Direction of Steps Probability 
x n one 
J} Jr FV 
Ky SrA * wu 
L- (n+ j, + kj) \7 ere 
y m one 
1) 5V AV an 
kK, uf/rA + y 
L- (m+ i, +k,) u/A tuty 
Zz R one 
1,5 v/u Fv 
J5 vV/i * Vv 
L- (£ +i, + j,) vV/rX tu ty 


where 1 = 1) +415, J) =], + Jog, kK =k, +k. Then 





P(Q) = Dalila eT a) Cpragras Ne ay 





cee (a i - (mt tk.) (yd aie 
At AtutY uty At 


bv, La (eti 35) > z 
ao) 22 ay? (Saag? 2 (4 2) 
~ —. Atuty, lL ,Atuty,y 
PAO) ee sap diy Kop 2 pp) rag ) (=a ) 














(Atutvyk | Vv yore U a rv ie 
A*tu Atuty AU AtutY 


where Q = Q(i,j,k,2,m,n) and r(i,j,k,£,m,n) is the number of paths 
from the origin to the point (L,L,L) which have i,),k,2,m,n points 


in common with boundaries A, B, C, D, E, F, respectively. 
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B. TRANSFORM OF THE DISTRIBUTION 


Fran equation (4-1) the following Laplace transform can be 
tabulated: 


Random Variable Laplace Transform 


A, (utv/utvts) i 

= teen 
(oaky/amnunnedhes 

D, (v/vts) * 

EM (u/uts)™ 

E. (A/\+s)™ 

G. (atutv/Atutves) 9 


The Laplace transform of the total replenishment time, T, is: 


f(s) =Ele™) = 2 Ble Io) p@ 





ie. - aie Apneee, Tate Assur. 3ilh 
foals) = =  xr(i,j,k,2,m,n) (Eby? mao 
1,J pie 


(Attys i Qtuivts) J (Atutvts) k (4-3) 
Ltvrs htt \+u+s 


ao ats R ,Atutves\M ,Atutvtes,n 
( | re it 
Vts ( Uts yf A+S 
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V. COUNTING TECHNIQUES 


A. DISCUSSION 

An expression for r(i,j,k,2,m,n) has not been developed in this 
paper. However, this chapter demonstrates several counting techniques 
for the two supply ship, two dimensional random walk which may be use- 
ful in the eventual determination of r(i,j,k,2,m,n). Additionally, 
for the three dimensional model an expression is found for the number 


of paths proceeding from the origin to (u,v,w) which touch no bound- 


i 


e = nl/r! (n-r) ! 


aries. In what follows, the binamial coefficient, ( 


B. TWO DIMENSIONAL COUNTING 
1. The Number of Paths Not Touching Boundaries 
First a method will be formulated for counting the number of 
paths proceeding from (0,0) to (u,v) which do not touch or cross either 
boundary y= x +a, (a > 0), or boundary y=x+b, (b < 0). This 
number is equal to the total (ae) paths going to (u,v) less the number 
of paths which touch or cross either boundary. By the Reflection 


Principle (see page 70 of [1)), the number touching or crossing y=x+a 


aod 
v-a 
in both of these quantities paths touching or crossing both boundaries 


is ) and the number touching or crossing y = x +b is Or)» Since 
are included, the number of such paths must be added back. This quantity 
may be computed by dividing it into two mutually exclusive categories -- 
those paths which touch or cross boundary y = x + a first, and those 
paths which touch or cross boundary y = x + b first. Such paths will 

be counted by a double application of the Reflection Principle. In 


Figure 3 consider path 1 which falls into the first of the two categories. 
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The first reflection is with respect to line y = x ta. Both the 
boundary y = x + b and that part of path 1 beyond its first cammon 

point with y = x + a are reflected. The reflected path, path l', 

ends in point (v-a, uta). The cammon points of path 1 and y=x+b 

are reflected into ccmmon points between path 1' and the line 

y = x + 2a-b (this being the reflection of line y=x +b). Therefore, 
the number of paths in the first category is equal to the number of 

paths ending in (v-a, uta) and touching or crossing the line y = x + 2a-b. 


The latter number is 


Vaatira Uv 


ia) é eas! (5-1) 


by a second application of the Reflection Principle. Similarly the 


number of paths going to (u,v) after touching y = x + b first is 


‘oats | ™ 


‘eam is Such a path is path 2 in Figure 3. Therefore the number 


of paths to (u,v) touching neither y=x+anory=xv+b is 


is +, _ /Utv uUtv utV , 
Oy) 7 Gate? 7 Soap! * Goes? 4 aie? (5-2) 


By extending the boundary reflection employed in the preceding 
paragraph it is possible to campute the number of paths proceeding to 
ayy) whieh touch or cross y = x + a, a > 0, first. y = x +), b << 0, 
second, and y =x +c, c > a, third, before going to (u,v). In Figure 
4 since y = x + a is the first boundary to be struck, reflect y=x+b 
and y = x + c with respect to y = x +a yielding y = x + 2a-b and 
y = x + 2a-c, respectively. Since y = x + 2a-b is the next boundary 
to be struck by the reflected path, reflect y = x + 2a-c with respect 


to y= x + 2a-b, yielding y = x+2a - 2b+c. The twice reflected path 


ley, 





| 
Me 
r (u-atb , vta-b) 
ah 

| 


Figure 3 
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(v-atb-c, uta-btc) 


2a-2b+c 


2a-b 





2a-c 


Figure 4 


ES 


will then strike y = x+2a - 2b+c and then proceed to (u-atb, vta-b). 
The number of such paths is, by the Reflection Principle, 


u-atb+v+a—b oh (atv 


ee RSD ee) ~ ‘weatb-c (5-3) 


This technique may be extended to multiple reflections of any number 
of similar boundaries. 
2. The Number of Paths Touching a Boundary d Times 
Consider now a technique for counting the number of paths to 
(u,v) which touch boundary y = x + a, (a > 0), d times, but do not 
cross it, and do not touch or cross boundary y=x+t+b, (b< 0). This 


number is equal to the number of paths proceeding to (u,v) which touch 


y=x+a d times without crossing it less the number of paths touching 


x +a d times without crossing it that do touch or cross y=x tb. 


y. 
In [4] Milch and Waggoner developed a new technique called the Telescope 


Principle for counting the number of paths proceeding to (u,v) which 


touch y +x +a d times and do not cross it. This number is equal to 


(utv—d (utv—d) 


a ~ * Uta (5-4) 


To campute the number of paths to (u,v) which touch y=x+ta 
d times, do not cross it, and touch or cross y = x + b, it is necessary 
to examine several disjoint subsets of that set of paths. 
If the number of paths in the set is N then 
N = N, + N., + N, where 


il Z 3 
the number of paths among these paths that touch 


ZA 
ll 


or cross y = x + b before touching y=xt+a 
N.. = the number of paths among them that touch y=xt+a 


d times without crossing it before touching or crossing 


y=xtb 


20 


(péq+n!q-A) 


(q4n/q-A) 


(a‘n) 


7 


G smMbTyZ 





ZL 


cell 
N, — (the number of paths among them that touch 
i=l 
y = x + a exactly 1 times and do not cross it before 
touching or crossing y = x + b) 


To find N, consider a path to (u,v) that crosses y=xt+b 


al 
first, then touches y=x+a da times without crossing it. Such 
a path is shown in Figure 5. If the path beyond its first point of 
contact with y = x + b is reflected about y = x + b, the reflected 
path (dotted line) will proceed to (v-b, utb) and touch the reflection 
of y=xt+a, i.e., the line y = x + 2b-a, d times without ever 
crossing it. The Telescope Principle is then applied to the reflected 
path. Therefore the number of such paths is equal to 

[number of paths to (v-b, u+b-d) hitting neither y = xta 

nor y = x+2b-a-d] 

~ [number of paths to (v-b, ut+b-d) hitting neither y = xta 
nor y = x+2b - a-d+l) 


This number is obtained from (5-2) yielding the formula for N 


L: 
m=, uty _ iv-dy Bee utv-d 
Na a aan tae UC 42a-b-1? 
utv—-d utv—-d utv—-d 
- oss IC -2at3b! 7 ac eere (5-5) 


Next the number of paths, N., proceeding to (u,v) which touch 


5! 
y=x+a d_ times but do not cross it then touch or cross y= x+b 
1s counted. Such a path is shown in Figure 6. If the Telescope 
Principle is applied to this path about boundary y = x + a, its 
equivalent path (dashed line) can be seen to meet the following 


Criteria: 


ZZ 


(A‘N) 





Q9 BaANbTy 





T-P+e 


Ze 


(1) 


(ii) 


it must touch but not cross line 
y=x+a+td-l first and then must 
touch or cross line y=x+t+betd. 
It may touch y= x + a+ d-l again 
after touching y=x+tb. 


Lt mast geeto (u=d7v)r 


The number of paths so constrained is 


where 


No 


ee ae —o Syyihriscl 
mA BSCR eDe ae ead) oe 
utv-d utv-d 
7 [C4 2a-b-1) 7 (yt2a-b? J (5-6) 


the number of paths to (u-d,v) which touch or 


cross y= x + a+ d-l first, then touch or cross 


ut+v-d 
vta-b-l’ * 


the number of paths to (u-d,v) which touch or 


y=x+b+d. From (5-1), A= ( 


cross y= x +a+d first,then touch y=x+be#d. 


utv-d) 
vta-b* *° 


the number of paths to (u-d,v) which touch or 


From (5-1) ,B= ( 


cross y = xX t+ a+ d-l first, touch or cross 


y = x + b +d second, then touch or cross 


utv—d 
v-2atb-dtl’ * 


the number of paths to (u-d,v) which touch or 


y=xtatd. From (5-3), C= ( 


cross y = x + a+td first, touch or cross 
y = X + btd second, then touch or cross 
y = X + atd again. D must be added since 


it has been subtracted in both B and C. 


utv—-d 


BE UE Dree Ne ead 
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Those paths which proceed to (u,v) and touch y= x +a (d-l) 
times before touching y = x +b, touchy =x+t+aée a total of d times, 
and do not touch y = x + a + 1 are counted next. Figure 7 shows a 
path of this type. Applying the Telescope Principle to this path 
about boundary y = x + a its equivalent path (dashed line) can be 
seen to be constrained as follows: 

(1) it must touch boundary y = x + a+ d-2 first, 

touch y = x + b + d-l second, and then touch 
Va= soerea tec. 

(ii) it mast go to (u-d,v). 

(iii) it may not touch y = x + a +d-l until it has 
Louewee y =x + b + d-l. 
(iv) it may not touch boundary y=x+t+atd. 


The number of such paths is equal to 


i. pen ee PV =a) 2 Paved 
Oe Sop Se a 
utv—-d utv—-d 


Es a eeb a (a+2a-b) | (5-7) 


where 
A' = the number of paths to (u-d,v) which must touch 
y = xtatd-2 first, touch y = x+btd-l second, then 
7 _ a SCI 
touch y = xtatd-1l. From (5-3), A' = (74 ai) 
B' = the number of paths to (u-d,v) which must touch 


y = xtatd-l first, touch y = x+btd-1l second, then 


(utv—d ) 
v-2atb-d+l1° ° 


C' = the number of paths to (u-d,v) which must touch 


touch y = xtatd-l. From (5-3), B' = 


y = xtatd-2 first, touch y = xtbtd-l second, then 


utv-d ) 


= es ( 
touch y = xtatd. From (5-3), C ( oatb-d+l 


Zo 


(A‘n) 





L sanbty 


q- T+P-q- 


Y 
‘e 


her 


C-P+e 
T-P+¥ 
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D' = the number of paths to (u-d,v) which must touch 
y=xtatd-1 first, touch y=x+tbt+d-1 second, then 
touch y=xtatd. D' must be added since it has 


been subtracted in both B' and C'. Fram (5-3), 


(atv-d ) 
v-2Zatb-d‘ ° 


When the numbers of those paths proceeding to (u,v) which 


D' = 


touch y=xta 1,2,...,(d-2) times, respectively, before touching 
y=xtb, touch y=xta a total of d times, and do not touch y=xtatl 
are computed similarly, each is found to be equal to (5-7). 


Therefore the total number of paths in N, is 


utv—-d utv—-d 


utv—-d 
SHRI SeAl! + (d-1) (e2a—b? 


ut2a-b-2’ a 


N, = (d-1) ( 


This may be written as 


2 ee . (neta 
- (6-1) (ab-1) > Groacb?! = 
Now N = N, +N, + Nz, and from (5-5), (5-6), and (5-8) 
N = (uiv—-d — (Btv—-d) (uw-d 4 Cis 
rao) uta-b u-2at3b U-2at3bt1 
Geman? ~ Gees) ~ ED (Cy zeny) BD a 
+ G1) (NDaeb-2) > Ga2acb-2)! 


Combining the results of (5-4) and (5-9), the number of paths 
to (u,v) which touch boundary y=xta d times, do not cross it, and do 


not touch boundary y=x+b is 


27) 


p(atv—-d) = (UtV—d) | _ pqatv-d a (Utv—d) | 


tard uta Draron uta-b 
ey cL utv—d utv-d utv—-d 
EG 2at3b) 7 “u-2at3be1)! 7 MGya-p-1) 7 rap?! Ga 
utv—-d vec _ utv—d ee uty cd 
ee LH +2a-b-1) 7 (i+2a-b’ | ce yt 2a-b=2) (a42a-b-1) 


C. THREE DIMENSIONAL COUNTING 

The number of paths from the origin to same point (u,v,w) of the 
three supply ship model which do not touch any of the boundaries of 
Chapter III will be counted next. This number is equal to the total 
number of paths fram the origin to (u,v,w), irrespective of boundaries, 
less the number of paths which touch at least one of the boundaries. 
The total number of paths from the origin to (u,v,w) irrespective of 
any boundaries is the number of possible combinations of u steps in 


the x direction, v steps in the y direction and w steps in the z 


direction. This number is the multinomial coefficient, ue x ) = 
tA 

utvtw, _ nee 

ere, (utvtw) | /alviw! 


In order to count the number of paths from (0,0,0) to (u,v,w) which 


contact boundary planes, it will be necessary to extend the Reflection 


Principle to three dimensions. By the Reflection Principle (a ; a7 
wta 7V¥,U-a 
os 
ileal ), and mae ) are the number of paths that have at least 
v-a ,uta,W u,w-a,vta 


one cammon point with the planes x=zta, y=xta, and z=yta, respectively, 


utvtw 


among the ( 
u,v ,W 


) paths connecting the origin and the point (u,v,w). 
If N' is the number of paths from the origin to (u,v,w) which touch 
at least one boundary plane then 
N'=E'+F'+G'-H'-I'-g'+k' 


where 
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E' = the number of paths touching x =zZz+M 
F' = the number of paths touching y=x+N 
G' = the number of paths touching z=y +P 
H' = the number of paths touching x=z+Mandy=x+t+N 


I' = the number of paths touching x=z2+Mandz=ytP 


x+Nandz=ytP 


J' = the number of paths touching y 
K' = the number of paths touching x =z+M, y= x +N, andz=ytP. 


By direct application of the Reflection Principle 


,  fahv-w 
F! = aig 
V-N ,UtN ,W 
a = (atv 
Uu ;w-P /VtP 


H' is the sum of two numbers; the number of paths which strike x =z+M 
first and the number striking y = x + N first. I' and J' are the sum of 
Similar terms. The formula for J' is explained below and H' and I' are 
listed in tabular form. The reflection of the plane y = x + N about the 
plane z=y+Pis x=z- (N+P). Any path proceeding to the point 
(u,v,w), initially touching the plane z = y + P (one or more times), and 
then proceeding to the plane y = x + N, is reflected beyond its first 
point of contact with z = y + P about that plane. The reflected path 
touches or crosses the reflection of y = x + N, i.e., the plane 


x = 2z- (N+P), and proceeds to the point [(u-N), (v-P), (Ww+N+P)]. 


utvtWw 


The number of such paths is eee 


The reflection of plane z = y + P about y=x+Nisx=2zZz- (N+P). 
Any path proceeding to the point (u,v,w), initially touching the plane 


y = x + N (ome or more times), and then proceeding to the plane z=y +P 


Zo 


is reflected beyond its first point of contact with y = x + N about 


that plane. The reflected path touches or crosses the reflection of 


Z=y +P, i.e., the plane x = z - (N + P), and proceeds to the point 


Hog = OMNES Be 


(v +N), (w+ P)]. Therefore J' 


i ay ) + (atvew 
u-N,v-P V+N ,wtP* * 


H' and I' are similarly computed and are tabulated below. 


Initial Plane 
Second Plane 
Reflected Plane 
Reflected Point 


Number of Paths 


utv-tw 
u-N,w-M 


a 


Initial Plane 
Second Plane 
Reflected Plane 
Reflected Point 


Number of Paths 


I' = ( 


ae 


| 


Table for H' 
Y=xXtN 
X+Z4+M 
z=y— (M+N) 


fu-N , V+M+N ,w-M] 


(ae 
u-N ,w-M 


UtVTW 
utM , Vv+N 


Table for I' 


Z=ytP 
x=zZ1M 
y=x- (M+P) 


[ut+M,v—(M+P) ,v+P] 


utvtw 


uM, w+P? 


( 


utv+w 
v—-P ,w-M 


X=Z+tM 
Y=xtN 
z=y— (M+N) 


[utM ,V+N ,w— (M+N] 


(ae 
utM ,VvtN 


y=x-— (M+P) 


[ut+M+P ,v—-P ,w—M] 


utV+Ww 
v—-P ,w-N 


( 


K may be computed in similar fashion by permiting all three boundary 


planes and applying the Reflection Principle twice after each permutation. 
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